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摘 要 : 针对 带 限制 的 开放 性 选 址 路 径 问 题 的 研究 , 考虑 模糊 需求 的 条 件 下 ,以 仓库 选 址 成 本 、 车 辆 行驶 距离 成 本 、 
机 会 损失 成 本 、 人 额外 距离 等 目标 之 和 最 小 化 的 要 求 下 建立 数学 模型 。 通 过 对 蘑菇 繁殖 算法 的 改造 ， 使 用 部 分 映射 交 
又 和 路 径 重 连 算法 代替 原 算法 中 父 代 更 新 方式 ; 在 邻 域 搜索 部 分 使 用 概率 法 进行 邻 域 选择 ; 使 用 随机 模拟 程序 对 设 
计 好 的 路 径 进 行 模拟 ， 计 算 因 服务 失败 而 产生 的 额外 行驶 距离 与 机 会 损失 成 本 。 在 保留 算法 原 有 特性 的 情况 下 ， 使 
其 成 功 应 用 于 组 合 优化 问题 。 最 后 ， 通 过 一 系列 算 例 测试 与 对 比 ， 验 证 了 模型 的 正确 性 与 有 效 性 以 及 混合 离散 蘑 茵 


繁殖 算法 的 计算 效率 和 优化 能 
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Hybrid discrete mushroom reproduction algorithm for solving open location-routing problem with 
fuzzy demands 


Liu Fan, Zhang Huizhen!, Zhou Xun 
(Business School, University of Shanghai for Science & Technology, Shanghai 200093, China) 


Abstract: This paper proposed the mathematical model which based on the open location routing problem with fuzzy demands 


considering with the sum of location cost, routing cost, lost opportunity and additional distances. For solving the mathematical 


model, this paper proposed a new swarm intelligence algorithm named hybrid discrete mushroom algorithm. The artificial 


wind replaced by partly mapping crossover and path relinking algorithm; And the local search strategy replaced by the method 


of selecting different neighborhoods with a certain probability; Stochastic simulation simulated the designed routes to compute 


additional distances and lost opportunity resulting from serving failure. The algorithm successfully applied to the 


combinatorial optimization problems under the condition of retaining the original characteristics of the algorithm. Finally, 


verify the validity of the mode and the performance of hybrid discrete mushroom reproduction algorithm through a series of 


tested instances and comparisons. 
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快速 、 精 确 的 车 辆 路 径 规划 与 合理 的 仓库 选 址 对 现代 物 


流 管理 有 着 重要 影响 ， 其 不 仅 可 以 节约 配送 成 本 ，T 


提高 客户 满意 度 及 企业 的 竞争 力 。 很 多 学 者 研究 指 


(facility location problem，FLP) 分 开 考 虑 时 ， 不 免 会 因为 所 考 型 ， 
虑 的 影响 因素 不 足 而 导致 管理 者 无 法 作出 最 优 决策 。 现 如 今 ， NSGA-II( 基 于 


两 者 的 综合 问题 选 址 路 径 问 题 (location-routing problem, 


LRP) 已 成 为 物流 管理 、 供 应 链 管理 及 运筹 优化 领域 的 有 


糊 限制 条 件 ， 


的 兴起 ,许多 企业 进行 战略 


liz 


户 需 求 、 服 务 时 间 等 的 不 确定 性 ,在 LRP 模型 中 加 入 各 种 模 
覆 成 带 模糊 需求 的 选 址 路 径 模型 中 ， 不 确定 性 
条 件 下 带 时 间 窗 的 选 址 路 径 模型 05， 多 元 联运 网 络 下 带 模糊 
需求 的 选 址 路 径 模 型 名 ， 应 急 物 流 系统 中 多 目标 模糊 选 址 路 
径 问 题 模型 站， 多 车 型 并 考虑 震 后 路 网 变化 情况 的 应 急 物 资 ”和 模拟 退火 算法 62 在 求解 LRP 问题 中 均 取 得 了 较为 满意 
动态 选 址 路 径 问题 模型 461。 此 外 ， 由 于 近年 来 第 三 方 物流 
组 ,为 了 提高 企业 核心 竞争 力 ， 
将 更 多 精力 专注 于 核心 业务 ， 而 将 物流 服务 外 包 给 第 三 方 公 


司 。 在 此 情况 下 ， 提 出 了 由 第 三 方 物流 配送 的 选 址 路 径 问 题 
一 一 开放 式 选 | 


j 且 可 以 件 ， 
: 将 车 中 对 带 模糊 限 秆 
辆 路 径 问 题 (vehicle routing problem，VRP) 与 设施 选 址 问题 中 ， 


究 热 点 。 ”配送 的 实例 。 
根据 实际 需求 , 在 标准 LRP 问题 中 加 入 不 同 的 约束 条 件 ， 
即 可 形成 不 同 的 LRP 扩展 问题 。 如: 根据 物流 配送 中 存在 客 。 生成 法 等 ) 仅 可 以 求解 小 规模 的 LRP 问题 ， 然 而 ， 从 实际 物 
流 配送 中 抽象 出 来 的 LRP 问题 往往 规模 较 大 , 精确 算法 难以 
间 内 给 出 问题 的 解决 方案 。 为 此 ， 在 有 限 


在 可 接受 的 


OLRP)。 虽 然 实际 配送 


引路 径 问 题 (Open location-routing problem, 


H, OLRP 也 会 面临 各 种 模糊 限制 条 


但 由 于 OLRP 是 近年 来 提出 的 新 间 题 ， 国 内 外 现 有 文献 


I 条 件 的 


建立 了 关于 OLRP 的 多 目标 混合 整数 规划 模 


OLRP 问题 研究 尚 少 ， 仅 在 文献 [15] 


利用 概率 进行 选择 更 加 安全 可 靠 的 路 径 进行 配送 ， 使 用 
基于 非 文 配 解 排序 的 遗传 算法 ) 和 NSDE( 基 于 非 文 


配 解 排序 的 差分 进 


结果 。 


LRP 属于 NP 难 问题 。 使 用 精确 算法 (如 分 支 定 界 法 、 列 


计算 时 


计算 时 间 内 能 给 出 近似 最 优 解 的 智能 优化 算法 已 成 为 LRP 
问题 的 主要 求解 方法 , 如 : 蚁 群 算 法 0 5 , SF EETEUSS T1771] 


本 文 将 模 糙 


需求 与 OLRP 相 结合 ， 提 出 了 带 模糊 需求 的 
开放 式 选 址 路 径 问 题 (Open location-routing problem with 


化 算法 ) 求 解 了 关于 汶川 地 震中 救援 物资 


ce 
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fuzzy demand, OLRP-FD), ifi 


EU 


E 够 很 好 解决 


凡 ， 等 : 


带 模糊 需求 的 开放 式 选 址 路 径 问题 的 混合 离散 蘑菇 繁殖 算法 


目前 第 三 方 物流 
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1 if r<dı 
配送 中 所 面临 的 客户 需求 的 不 确定 性 问题 ， 对 节约 企业 配送 TN nM 
CK RU SOERESGL. wih KERR RECTA s Ud p a] crd | 0 Ea A 
入 新 的 邻 域 搜索 方法 ， 设 计 求解 OLRP-FD 问题 的 离散 蘑菇 QE. died: 
繁殖 算法 ， 这 不 仅 为 OLRP 问题 的 求解 提供 了 一 种 新 的 解决 Ap 
方案 ， 而 且 拓宽 了 蘑菇 繁殖 算法 的 应 用 领域 。 算 法 中 的 邻 域 Du 
搜索 在 评估 搜索 父 代 蘑菇 邻近 区 域 的 基础 上 ， 生 成 子 代 孢子 ”1.2 问题 描述 
来 细 化 这 些 区 域 以 搜索 更 优 的 解 ， 数 值 实验 表明 该 算法 在 有 本 文 在 可 信 度 理论 的 基础 上 ， 构 建 带 模糊 需求 的 开放 式 
限 的 计算 时 间 内 可 为 OLRP-FD 提供 较 优 的 解 。 选 址 路 径 问 题 的 数学 模型 ， 由 选中 的 仓库 派出 车 辆 为 客户 提 
ut 供 配 送 服务 , 并 假设 : a) 每 辆 车 只 能 服务 最 多 一 条 线路 ; b) 仓 
库 和 车 辆 有 固定 的 容量 限制 ; (3) 配 送 车 辆 有 最 大 行驶 距离 限 
OLRP-FD 相 较 于 标准 OLRP 的 不 同 之 处 在 于 客户 的 需 W: c) 车 辆 服务 完 最 后 一 个 客户 后 , 并 不 会 返回 仓库 , 而 是 返 
求 是 模糊 的 ， 而 非 确定 的 ， 只 有 配送 车 辆 在 抵达 客户 时 ， 才 ” 回 第 三 方 物流 企业 车 辆 集散 中 心 ，d) 由 于 每 个 客户 需求 的 模 
能 知道 客户 的 具体 需求 。 因此 , 在 OLRP-FD 中 ， 当 配送 车 辆 。 糊 性 ， 可 能 会 存在 因为 车 辆 数 不 足 而 有 客户 无 法 被 服务 的 情 
在 抵达 客户 时 ， 可 能 会 发 现 因为 车 辆 的 剩余 容量 不 足 而 无 法 。 况 ， 此 时 会 产生 相应 的 机 会 损失 成 本 。 
服务 当前 客户 ， 此 时 ， 车 辆 需要 返回 仓库 重新 装载 ， 之 后 再 于 客户 需求 是 模糊 的 这 一 情况 ， 文 章 利用 可 信和 度 理论 
可 到 未 被 服务 的 客户 进行 服务 。 OLRP-FD 中 模糊 的 客户 需求 ”来 判断 当前 车 辆 以 及 仓库 是 否 能 服务 下 一 个 新 的 客户 。 假 设 
不 仅 会 导致 配送 车 辆 返回 仓库 重新 装载 带 来 额外 的 行驶 距离 ， ”车辆 与 仓库 的 容量 分 别 为 0 和 P， 当 车 辆 或 仓库 服务 完 计划 
而 且 会 因 配 送 车 辆 预先 配备 不 足 而 无 法 服务 部 分 客户 ,产生 ”线路 中 的 第 e 个 客户 时 ， 当 前 车 辆 的 剩余 容量 为 
一 定 的 机 会 损失 成 本 。 e : i em 
本 文 使 用 可 信 度 测量 理论 来 构建 OLRP_FD 的 数学 模型 ， 070 200 2450 2,00 n (mum, 当前 仓库 的 剩余 
模型 在 满足 车 辆 容量 、 最 大 行驶 距离 等 条 件 下 ， 要 求 不 仅 最 oun. -—' 
AouR AME EMHE dE a te 0 0 T 
离 成 本 、 未 服务 客户 的 机 会 损失 成 本 ), 而 且 最 小 化 车 辆 由 于 ”一 个 客户 需求 不 超过 车 辆 容量 和 仓库 容量 的 可 信和 度 可 分 别 表 
服务 “失败 ”而 额外 行驶 的 总 距离 。 示 为 式 (4) 和 (5): 
11 可 信和 度 理 论 Cr(d.-i€ Qo}=Cr{(G, «i quads e+1-q2ed3, eri~gue) <0} (4) 
模糊 集 理论 最 先 由 Zadeh20 通 过 隶属 度 函 数 提 出 , 之 后 ， Cr(d. € Pe} = Cr((di pod e+1— pze ds. eri-pre) <0} (5) 
由 Kafumann??!fll Zadeh22 完 善 了 测量 模糊 事件 的 模糊 变量 进而 ， 由 式 (3) 可 将 crtd scj 和 Crid P 分别 表示 为 
可 能 性 测量 理论 ， 但 其 缺乏 自 对 偶 性 。Liug 为 定义 其 自 对 偶 。 下 式 (6) 和 (7): 
性 ， 提 出 了 关于 模糊 变量 的 可 信 度 测量 理论 。 ro WT 
假设 ( €. P(O), Pos ) 为 可 能 性 空间 , A 是 PO 中 的 一 个 事件 ， 5 uoti. 
Pos{ 和 表示 事件 4 发 生 的 可 能 性 ， 那 么 4 的 必要 性 (也 就 是 ge adie i 
事件 A 发 生 的 必要 性 ) 则 为 Nec{4}=1- Pos{4*} 。 记 4 为 事件 4 duds drertdaer gd "PD cus 
的 补 集 , 进而 可 得 事件 4 的 可 信和 度 Cr{A}=1/2(Nec{A}+ Post A53) 。 0 
本 文中 将 客户 需求 使 用 三 角 模 糊 变 量 进行 表示 为 di op a SA 
d-(dud»ds) , $F, di 和 4d; 分别 是 客户 需求 的 上 下 界 ， 客 3es1É qe 
户 的 需求 一 定 在 这 个 范围 之 内 ;qd 表示 隶属 度 函数 值 为 1 的 oops 
时 候 所 对 应 的 客户 需求 。 这 三 个 参数 一 般 由 相关 数据 或 者 经 
验 主观 估计 得 到 。 假设 一 个 客户 的 实际 需求 为 , 隶属 度 函 数 T "mdp Mn 
如 图 1 所 示 , 其 可 能 性 , 必要 性 , 可 信和 度 可 分 别 由 公式 (1)-(3) Dr (pre dien tdzen pd) PP CT) 
计算 所 得 。 dS pu Ue DNO. oe pe di pie 
2*(p2ze-d2e+1+d3,e+1— pie) 
mo) 1 daer1< pie 
Cr. 为 评估 配送 车 辆 或 仓库 是 否 为 某 个 客户 提供 服务 的 
参数 ， 当 Cr 为 0 时 ， 这 个 客户 必定 不 会 被 当前 车 辆 或 仓库 
服务 ; 当 Cr 为 1 时 , 该 客户 肯定 会 被 当前 车 辆 或 仓库 服务 。 
文章 设置 了 两 个 阔 值 DPI 和 A4PI1( DPI e[0,1} API e[0,1] ) 分 别 与 
udi 车 辆 和 仓库 的 可 信 度 进行 比较 , 当 Cr 大 于 所 设置 的 阔 值 时 ， 
图 1 模糊 变量 隶属 度 函 数 客户 会 被 当前 车 辆 或 仓库 服务 ;而 小 于 所 设置 的 阔 值 时 ， 则 
Fig. 1 Membership function of triangular fuzzy variable 不 被 当前 车 辆 或 仓库 服务 。 显 然 ，DPI 和 API 的 取 值 大 小 对 
i ff rdi 车 辆 计划 行驶 距离 和 额外 行驶 距离 有 着 重要 影响 。 当 DPI 取 
Po o acu d 值 较 低 时 ， 车 辆 的 容量 会 被 尽 可 能 使 用 ， 此 时 ， 路 线 中 靠 后 
di-di ` 的 客户 就 极 有 可 能 服务 失败 ， 从 而 导致 车 辆 要 返回 仓库 重新 
9 Ue 装载 ， 产 生 额 外 的 行驶 距离 ， 当 DPI 取 值 较 高 时 ， 意 味 着 仓 
i jf < 库 派出 的 车 辆 数 较 多 ， 那 么 计划 行驶 距离 就 会 增加 ， 而 服务 
E E D 失败 的 客户 就 会 减少 ， 产 生 的 额外 距离 也 随 之 减少 。 当 API 
dad ` 取 值 较 低 时 ， 表 明 仓库 会 尽 可 能 的 消耗 自己 的 容量 去 服务 客 
grew 户 ， 不 仅 有 可 能 导致 其 容量 溢出 ， 而 且 有 可 能 导致 仓库 预 留 
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di 第 7 NAPR, JeJ; 
Qk: 第 上 辆 车 的 容量 ，keK ; 


Pi: 第 i 


个 仓库 的 容量 ，iel ; 


O: 第 i 个 仓库 的 开放 成 本 ，iel ; 
Fa 第 大 辆 车 的 启用 成 本 ， keK; 


ADi 第 上 辆 车 行驶 的 额外 距离 ， KeK ; 


Ci: 点 


i 到 点 j 的 单位 距离 行驶 费用 ，ieV, jev; 
TDi: 车 辆 的 最 大 行驶 距离 ，keK; 

Dij: 点 1 到 点 j 的 欧 氏 距离 ，ieV, jeV; 
Bi 未 服务 客户 7 时 的 损失 成 本 ， jeJ7 。 


决策 变量 设置 如 下 : 


1 
Xi — 
0 


车 辆 k 从 i 点 行驶 到 j 点 ,VieV,vijeV,vkeK 


其 他 
y -f 仓库 i 服务 客户 Viel, YjeJ 
lo 其 他 
j= f 仓库 i 开放 ,vie 了 
0 其 他 


在 上 述 假设 条 件 及 模型 参数 和 决策 变量 设置 的 基础 上 ， 
OLRP-FD 的 数学 模型 可 构建 如 下 : 


min >， OiZi+ YX X FiXi + YXYX CiDiXii 
icl icl jc] keK icV je] kcK 


min) ADS (0 - Y Y Xi)B; 
keK jel ieV keK 


Cr( 2 > dXasQO2DPI VkeK 
ieV jeJ 


CrC2 diYy < PZ)zAPI Viel 


jel 


X X< vjeJ 
ieV keK 


Wr—WritNXx <N-1  VijeJ;:VkeK 


> > Xa<l vkek 


iel jeJ 


> Xo) Xx=0 vjed;Vkek 
ieV ieS 


Ip Y XssieYs viernvjeJ;vVkek 


levy) 


X %=1 vjeJ 
iel 


> $ DiX ADi < TD, Yk € K 
V J 


ieV je. 


YXa-Zs0 viel,vkek 
J 


je 
YX2-X&s0 Viel 
jeJ keK 


Xie{0,1} VieV;VjeS;VkeK 


(8) 


(9) 


(10) 


a1) 


(12) 
(13) 
(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


Q1) 
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的 配送 车 辆 不 足 ， 部 分 被 安排 给 此 仓库 服务 的 客户 无 法 被 提 Yie(01] ViervjeJ (22) 
供 服 务 , 从 而 产生 较 大 的 机 会 损失 成 本 ; 当 4P7 取 值 较 高 时 ， Ze(01) Viel Q3) 
仓库 的 容量 对 客户 的 需求 安排 更 加 谨慎 ， 因 而 出 现 容量 溢出 Wre{OUN} ieJ;keK (24) 
的 情况 比较 少 。 本 文通 过 随机 模拟 程序 求解 相应 的 额外 行驶 其 中 ， 目 标 函 数 (8) 表 示 最 小 化 总 成 本 ,包括 仓库 和 和 车辆 的 固 
距离 ， 并 且 会 对 DPI 进行 灵敏 度 分 析 ， 求 得 最 优 参 数 设置 。 ” 定 启 用 成 本 、 车 辆 计划 行驶 总 距离 成 本 ; 目标 函数 (9) 表 示 最 
13 ”模型 小 化 由 服务 “失败 ”而 导致 的 车 辆 行驶 的 额外 距离 之 和 以 及 
模型 主要 参数 设置 如 下 : 未 被 服务 客户 的 机 会 损失 成 本 ; 模糊 机 会 约束 (10) 保 证 配送 
Lh CERA; 车 辆 对 所 有 客户 进行 访问 时 ， 车 辆 容量 与 客户 需求 满足 所 设 
J: 客户 集合 ; 置 的 可 信和 度 水 平 ， 模 糊 机 会 约束 (11) 表 示 仓 库 为 客户 提供 服 
JK: 第 三 方 物流 公司 车 辆 集散 中 心 点 ; 务 时 , 其 容量 与 客户 需求 满足 所 设置 的 可 信 度 水 平 ; 约束 (12) 
V: IDJ; 表示 一 个 客户 最 多 被 一 辆 车 服务 ; 约束 (13) 用 于 消除 子路 径 ; 
S: JKUJ; 约束 (14) 表 示 任 意 车 辆 最 多 被 启用 一 次 ; 约束 (15) 保 证 路 径 的 
Ki 车 辆 集合 ; 连续 性 ， 进 出 每 个 点 的 流量 相同 ; 约束 (16) 和 (17) 表 示 客 户 被 


某 个 仓库 服务 时 , 必定 有 路 线 将 两 者 相连 ; 约束 (18) 表 示 车 辆 
的 最 大 行驶 距离 限制 ; 约束 (19) 和 (20) 表 示 两 个 决策 变量 之 间 
的 关系 ; 约束 (21-23) 表 示 3 组 不 同 的 0-1 决策 变量 ; 约束 (24) 
为 消除 子路 径 所 用 的 辅助 变量 。 


2 ”求解 OLRP-FD 的 混合 离散 蘑 妆 繁殖 算法 


蘑菇 繁殖 算法 04(mnushroom reproduction optimization, 
MRO) 是 基于 蘑菇 生长 繁殖 特性 而 提出 的 一 种 新 型 群 智能 优 
化 算法 ， 已 在 连续 Mu ne Tue: 为 
了 进一步 拓宽 MRO 算法 的 应 用 领域 ， 本 文 将 其 与 贪心 聚 类 
算法 、 路 径 重 连 算法 mr 结合 ， 设 计 求 解 OLRP- 
FD 的 混合 离散 藤 冶 繁殖 算法 (hybrid discrete mushroom 
reproduction optimization, HDMRO). 

HDMRO 由 三 部 分 组 成 : a) 初 始 化 : 使 用 贪心 聚 类 算法 
生成 一 组 高 质量 的 初始 解 ， 以 便 算法 能 够 更 好 地 寻 优 迭代 ; 
b) 离 散 蘑菇 繁殖 算法 : 在 基本 MRO 算法 的 基础 上 ， 结 合 
OLRP-FD 问题 的 具体 特征 , 与 部 分 映射 交叉 和 路 径 重 连 算法 ， 
以 及 邻 域 搜 索 相 融合 ， 逐 步 使 问题 的 解 得 到 优化 ;c) 随 机 模 
拟 程序 : 通 过 模拟 客户 的 实际 需求 ， 计算 车 辆 服务 时 产生 的 
额外 距离 。 

2.1 基本 蘑菇 繁殖 算法 

蘑菇 繁殖 算法 09 通 过 模拟 蘑菇 繁殖 和 生长 的 机 制 ， 父 代 
蘑 要 通过 风 的 力量 来 传播 自己 的 孢子 ， 让 其 探索 各 个 繁殖 区 
域 ， 并 细 化 搜索 空间 ， 找 到 更 优秀 的 繁殖 区 域 来 进行 生长 繁 
殖 搜索 活动 。 基 本 芒 菇 繁殖 算法 最 初 是 用 来 解决 连续 优化 问 
题 的 ， 其 关键 步骤 为 : 在 进行 局 部 搜索 之 前 ， 根 据 每 个 菌落 
的 平均 适应 度 值 (4ve(i)) 和 所 有 菌落 的 平均 适应 度 值 Tags F 
断 Avg) Tos! c 与 Tag 的 大 小 关系 (c 为 一 个 固定 的 阔 值 ， 
c e[L101), 决定 是 否 需要 删除 一 些 菌 落 。 若是 前 者 大 于 后 者 ， 
则 不 进行 父 代 莽 菇 的 删除 , 各 小于, 则 对 父 代 蘑 菇 进行 更 新 ， 
选取 更 加 优秀 的 父 代 细 化 搜索 空间 , 以 便于 搜索 到 优秀 的 解 。 
对 父 代 蘑菇 进行 更 新 时 ， 通 过 构造 人 工 风 来 进行 钨 子 之 间 信 
息 交 换 更 新 。 第 7 个子 代 戎 菇 人 工 风 的 初始 计算 公式 为 (25): 


Mu ) 


TA 


Mov”! 2(X? — Xj)x( )-? x rand (0,6) x rs + rand(—s, s) (25) 


其 中 , XR XCAY2J7958 i iE ZAR; B (Bello) ) 
为 风 的 密度 ，; 5(5e(0,D ) 为 调整 风 的 方向 参数 ;rs (7se(0,D ) 
为 控制 风 的 步 长 的 参数 ; s 为 搜索 半径 , 其 通常 根据 具体 问题 
定义 。 


除 上 述 步骤 外 ，MRO 算法 中 通常 采用 下 述 式 (26) 进 行 邻 
域 搜索 : 


X, = Xi &rand(-s,s) (26) 
Ep, Xy HE iA ACROSS AE ETT SR JT RC 
基本 MRO 算法 主要 通过 式 (23) 和 (26) 进 行 迭 代 搜 索 ， 有 具 
体 迭 代 优化 步骤 参见 文献 [14] 。 
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2.2 初始 化 


2.2.1 解 的 表示 


本 文采 用 
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然 数 编码 。 假 设 m 个 仓库 的 编号 分 别 为 1， 


25 ess m n 个 客户 的 编号 分 别 为 m+ , m+2, ..., MHNA, 
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类 客户 的 顺序 ， 随 机 构造 一 条 配送 路 线 。 
2.8 改进 蘑菇 繁殖 算法 
2.3.1 部 分 映射 交叉 和 路 径 重 连 算法 

为 了 将 蘑菇 繁殖 算法 很 好 地 应 用 于 LRP 问题 的 求解 , 提 


则 3 个 仓库 ，15 个 客户 的 OLRP-FD 的 算 例 的 可 行 解 可 编码 


为 图 2 所 示 的 形式 。 若 某 一 仓库 编号 之 后 无 客户 编号 ， 


则 表 


高 其 求解 性 能 , 本 文 利用 部 分 映射 交叉 P9 和 路 径 重 连 算法 PC 
蔡 换 基 本 蘑菇 繁殖 算法 中 对 父 代 蘑菇 利用 风 的 力量 进行 全 局 


示 该 仓库 不 开放 若 某 一 仓库 编号 之 后 为 客户 编号 ， 则 表示 


该 仓库 开放 ， 而 
配送 路 线 数 。 一 辑 


车 从 对 应 编号 的 仓库 4 


编号 为 仓 请 


对 相应 客户 进行 
配送 任务 。 妇 


8 现 的 次 数 即 为 从 该 仓库 出 发 的 
上 发 ， 从 左 到 右 依次 
FE 时， 该 车 辆 完成 


J2 所 示 ， 仓 库 1 关闭， 仓库 2 和 3 开放 ; 从 


仓库 2 出 发 的 2 条 配送 路 线 分 别 为 2 一 5 一 7 一 4 一 10 和 2 一 6 


一 12 一 14 一 1 
一 9 一 19 一 10 


3 一 18; 从 仓库 3 H 


出 仓库 所 要 月 


一 8。 这 样 的 编码 方式 简单 直观 ， 能 很 好 地 解读 
户 ， 每 辆 车 所 行驶 的 路 线 。 


2.2.2 初始 化 


初始 可 行 解 的 优 劣 不 
法 的 收敛 速度 。 为 了 得 至 
聚 类 算法 进行 初始 化 ， 

1) 聚 类 客户 

利用 贪心 思想 ， 根 据 车 软 
及 客户 之 间 的 距离 对 客 


a) 生成 一 个 空 聚 


响 解 的 质量 ， 而 且 
| 较 优 的 初始 可 行 解 ， 本 文 使 


类 ( 聚 类 的 容量 等 


Fig.2 The representation of solution 


发 的 配送 路 线 为 3 一 17 一 16 


搜索 的 处 理 方式 (23)。 有 具体 方法 为 : 首先 , 将 第 i 个 菌落 的 父 
代 蘑 茹 与 全 局 最 优 解 或 第 上 个 菌落 的 父 代 蘑菇 进行 部 分 映射 
交叉 , 得 到 新 解 ; 然后 使 用 路 径 重 连 算法 对 当前 迭代 菌落 ;的 
父 代 与 部 分 映射 交叉 所 产生 的 解 进行 连接 ， 搜 索 两 者 之 间 更 
优 的 解 。 

部 分 映射 交叉 : 产生 随机 数 re[0,1] , 若 r<05， 选择 全 局 
最 优 解 与 第 i 个 菌落 的 父 代 藤 阁 进行 交叉 ; 若 r>0.5 , 选择 第 


会 影响 算 


1 信心 


户 的 模糊 需求 ， 以 
户 聚 类 。 一 次 聚 类 过 程 可 简要 描述 为 
F 车 辆 的 容量 )， 并 从 未 被 


聚 类 的 客户 中 随机 挑选 一 个 客户 加 入 该 聚 类 ; b) 从 未 被 聚 类 


的 客户 中 ， 选 择 距 新 力 


"最 近 的 客户 ， 并 计算 


此 客户 的 模糊 需求 与 车 辆 剩余 容量 的 
给 定 的 DPI， 则 将 该 客户 加 入 当前 聚 类 ; E 
考虑 此 客户 ; c) 8 
聚 类 的 客户 最 近 的 
直到 当前 聚 类 的 匀 


可 信和 度 ， 若 其 值 大 于 所 
则 ， 当 前 聚 类 不 
BK b) ， 继 续 选 择 离 最 后 一 个 被 加 入 
并 判断 其 是 否 可 以 加 入 当前 
不 足以 服务 其 他 新 客 


is 
m, 


户 时 ， 一 次 聚 


类 完成 。 重 复 步 又 a)j~c) 继 续 生成 新 的 聚 类 ， 直 到 所 有 客户 都 


limi 


被 聚 类 为 止 。 


DE 


iu He 4 NU PR 
首先 , 根据 每 个 聚 类 a 
ib. 然后， 由 式 (28) 和 (29) 分 别 计算 每 个 仓 


式 (27) 计 算 聚 类 重 
库 到 客户 聚 类 重心 


的 距离 之 和 U, 和 排序 指标 W。 最 后 , 按 W, 


值 从 大 到 小 8 


的 顺序 


对 仓库 进行 排序 ， 开放 W 值 最 大 的 仓库 ， 
户 聚 类 分 配给 它 ， 接 着 分 本 


仓库 剩余 容量 不 足以 月 


序列 中 下 一 个 仓 | 
某 一 客户 聚 类 分 本 


了 聚 类 都 被 分 配给 仓库 为 
给 某 一 仓库 时 ， 需 要 考虑 该 聚 类 的 总 需求 


将 离 它 最 近 的 客 
离 它 次 近 的 聚 类 ， 依 此 类 推 
R 务 任意 未 被 分 配 的 聚 类 时 ， 开 启 排序 


NA 
| 


EF。 将 


量 与 仓库 剩余 容量 之 间 的 可 信和 度 值 是 否 大 于 给 定 的 API 值 ， 


若 大 于 ， 则 将 聚 类 
的 聚 类 是 否 满足 该 条 件 ， 直 到 仓库 


为 止 。 


3) 构建 配送 线路 
对 每 一 客户 聚 类 及 为 其 分 配 的 仓库 ， 按 照 先 仓库 ， 


分 配给 该 仓库 ， 否则， 检验 下 一 未 


<“ )，c 为 聚 类 
m 


Ui Z3 (X-X) - (i -Y)?.,i el 
ceJ 


D^) RC 
不 足以 服务 任意 聚 类 


个 菌落 的 父 代 蘑菇 与 第 i 个 菌落 的 父 代 蘑菇 进行 交叉。 实验 
证 明 本 文 所 采用 的 部 分 映射 交叉 算法 可 以 有 效 加 快 算法 收敛 
速度 。 
图 3 给 出 仓库 数 和 客户 数 分 别 为 2 和 6 的 OLRP-FD 实 
例 的 部 分 映射 交叉 过 程 。 如 图 3 所 示 : 首先 ， 随 机 选择 父 代 
P1 和 P2 中 的 两 个 位 置 4 和 7， 将 父 代 个 体 第 4 到 第 7 位 中 
所 有 信息 复制 到 其 相应 的 子 代 个 体 Childl 和 Child2 中 ( 子 代 
个 体 Childl 和 Child2 中 其 余 位 置 用 * 表 示 )， 得 到 Child1(*) 
和 Child2(*)， 以 及 一 组 映射 (4 一 9，1 一 6，5 一 2，7 一 4); B 
次 ， 从 左 到 右 依次 扫描 个 体 P1 和 P2， 若 父 代 P1(P2) 中 第 i 
个 位 置 上 的 编码 未 在 子 代 Child2(Child1) 中 , 则 将 其 复制 于 子 
代 Child2 的 相应 位 置 上 , 得 子 代 Child1”(*) 和 Child? (*); 再 
次 ， 利 用 所 得 映射 填充 子 代 ChilaU (*) I Child2”(*) 中 的 剩余 
空位 ， 如 : FA Child1”(*) 中 的 空缺 位 置 1，2 和 3 本 应 填 入 
父 代 P2 中 相应 位 置 的 编码 1，7 和 5, 但 由 于 Child H E 
存在 这 些 编码 , 则 根据 映射 规则 (1 一 6, 7 一 4 一 9, 5 一 2) 将 1, 
7 和 5 分 别 调整 为 6,9 和 2, 得 ChildU , 同 理 , 可 得 Child2’; 
最 后 ， 检 查 Childi’ $ Child?’ $ 1 个 位 置 的 编码 ， 若 Childl 
和 Child2'28 1 个 位 置 的 编码 为 客户 编码 (如 图 3 所 示 , Child1 
和 Child? $ 1 个 位 置 的 编码 分 别 为 6 和 5) , 则 随机 选取 任意 
一 个 仓库 (如 图 3 所 示 , 选取 仓库 1) 与 其 交换 , 得 最 终 解 Childl 
和 Child2 。 


Joar 
e| 


图 3 部 分 映射 交叉 
Fig.3 Partly mapping crossover 
路 径 重 连 算法 : 路 径 重 连 算法 主要 是 对 两 个 解 (作为 当前 
种 群 父 代 孢子 的 当前 解 X 和 部 分 映射 交叉 得 到 的 目标 解 XT) 
言 息 的 交换 。 算 法 具体 过 程 可 被 描述 如 下 : 
a) 计算 初始 解 玉 与 引导 解 邓 之 间 的 汉 明 距离 HD, 并 令 


Q7) 


Q8) 


Q9) 


后 聚 


标 解 并 = 和 
b) zi HD>0 H. HD 为 奇数 ， 依 次 从 左 到 右 比较 X E; x 
相同 位 置 的 编码 , 并 对 不 同 的 编码 进行 交换 ; 关 HD>0 H. HD 
为 偶数 , 依次 从 右 到 左 比 较 凶 与 承 相 同位 置 的 编码 ， 并 对 不 
同 的 编码 进行 交换 ; 若 HD=0， 则 输出 目标 解 Y。 交 换 过 程 
可 概述 为 : B X 5 X E a 个 位 置 的 编码 不 同 ， 在 X PER 
到 与 系 (g 的 编码 相同 的 位 置 b, S XXa), X(a)-X'(a); 

c) 更 新 HD: = HD-1， 并 交换 当前 解 邓 与 引导 解 X 
(X'eoX), GTP b). 

路 径 重 连 算 法 可 以 免 去 大 多 数 连续 优化 的 算法 在 求解 离 
散 问题 时 ， 解 的 信息 从 离散 到 连续 ， 再 从 连续 到 离散 的 转换 
过 程 ， 有 效 防止 了 信息 的 丢失 。 
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2.3.2 邻 域 搜索 


刘 JU, €i 带 模糊 需求 的 开放 式 选 址 路 径 问 题 的 混合 离散 蘑菇 繁殖 算法 


为 了 进一步 提高 蘑菇 繁殖 算法 在 离散 优化 问题 中 的 寻 优 


能 力 , 本 文 在 基本 蘑菇 繁殖 算法 的 
题 的 特点 ,选择 了 单 点 交换 、 扣 


3 种 邻 域 以 相同 概率 对 解 进 行 搜 索 。 


[ 交换 位 置 


基础 上 


图 4 单 点 交换 


Fig.4 Exchange single point 


单 点 交换 : 随机 选择 两 个 位 


码 进行 交换 。 如 图 


17， 将 其 


解 的 可 行 性 ， 只 能 


, 结合 OLRP-FD 问 
入 和 逆转 3 种 邻 域 搜索 机 人 


e 


LURLJ ,将 其 相应 位 置 的 编 


4 所 示 ， 位 置 7 和 14 的 编码 分 别 为 6 和 


交换 后 ， 得 到 新 的 个 体 。 需 要 注意 的 是 : 为 了 保证 


进行 客户 与 客户 之 间 的 交换 和 仓库 与 仓库 


之 间 的 交换 ， 而 客户 与 仓库 之 间 不 能 进行 交换 。 
随机 选择 两 个 位 置 i 和 j，, 将 位 置 i 的 编码 插 
5 所 示 , 将 位 置 7 的 编码 17 插 


单 点 插入 : 


入 到 位 置 J 的 编码 之 前 。 如 图 
入 到 位 置 14 的 编码 6 之前， 得 到 新 的 个 体 。 
逆转 ; 随机 选择 两 个 位 置 i 和 


7 ,将 位 置 i 和 j 之 间 ( 包 括 


位 置 i 和 jj) 的 所 


14 之 间 的 编码 逆转 排序 后 ， 得 到 革 


编码 反 向 排序 。 如 


图 5 单 点 插入 


Fig.5 Insert single point 
[UN 


! 


图 6 逆转 


Fig.6 Inverse 


2.4 ”随机 模拟 程序 


如 前 文 所 述 ， 


由 于 在 OLRP-FD 


都 是 


图 6 所 示 , 将 位 置 11 和 
的 个 体 。 


模型 中 每 个 顾客 的 需求 


角 模 糊 变 量 表示 ， 所 以 ， 当 HDMRO 优化 好 仓库 选 


址 位 置 和 车 辆 行驶 路 线 之 后 ， 本 文 使 用 随机 模拟 0 的 方法 模 


拟 每 个 客户 的 实际 需求 , 进而 计算 机 会 损失 成 本 (由 于 仓库 
时 产生 的 机 会 损失 成 本 ) 和 车 辆 行 
j 按 照 规划 路 线 行驶 至 某 个 客户 时 , 车 
勿 量 无 法 满足 当前 服务 客户 的 需求 ， 从 而 车 辆 需要 
新 装载 货物 来 有 
离 ， 距 离 为 往返 当前 客 


辆 数 不 足 以 服务 所 有 客户 
驶 的 额外 距离 (由 于 车 软 


MAURE 
返回 仓库 


IT qu 


由 于 车 辆 在 抵达 当前 客户 之 


民 务 当前 客 


后 ， 发 


P, H 
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c) 重复 上 述 两 个 步骤 400 次 , 取 平 均 额 外 距离 和 平均 损 
失 成 本 为 最 终 的 额外 距离 与 损失 成 本 。 
2.5 HDMRO 算法 流程 

HDMRO 算法 的 具体 步骤 如 下 : 

a) 初始 化 参数 。 父 代 蘑 菇 数量 (种 群 数量 M)， 每 个 父 代 
产生 孢子 个 数 N， 算 法 最 大 办 代 次 数 G(g = 1D)， 靖 值 c。 

b) 调用 贪心 聚 类 算法 生成 初始 可 行 解 。 每 个 父 代 蘑 菇 通 
过 邻 域 搜索 生成 W 个 孢子 ， 计 算 每 个 个 体 的 适应 度 值 ， 选 取 
每 个 父 代 葡 问 和 其 产生 的 孢子 中 的 最 优 适 应 度 值 的 个 体 作为 
新 的 父 代 ， 更 新 每 个 父 代 蘑 茹 。 

c) 正式 开始 迭代 。 若 hvgO+Tes/c>7Tos ， 转 步骤 由 ; T 
N, FEIR e) 。 

d) 调用 部 分 映射 交叉 生成 一 个 路 径 重 连 算法 所 需 的 引 
导 解 ， 再 调用 路 径 重 连 算法 连接 当前 迭代 次 数 的 第 i 个 种 群 
的 父 代 解 与 部 分 映射 交叉 生成 的 引导 解 ， 搜 索 两 个 解 之 间 的 
区 域 。 更 新 每 个 菌落 的 父 代 蘑菇 和 全 局 最 优 值 。 

e) 邻 域 搜索 。 利 用 概率 法 对 三 种 邻 域 进行 选择 ， 父 代 蘑 
菇 根据 选择 出 的 邻 域 进行 繁殖 ， 产 生 孢 子 ， 搜 索 解 空间 。 更 
新 每 个 菌落 的 父 代 蘑菇 和 全 局 最 优 值 。 

Dg-gtl £i g«G, ER o, PERR: T9). fa 
出 当前 最 优 解 ， 并 转 步 又 g) 。 
g) 将 优化 好 的 路 线 使 用 随机 模拟 程序 进行 模拟 , 求解 出 
相应 的 额外 行驶 距离 和 未 被 服务 客户 带 来 的 机 会 损失 成 本 。 


初始 化 参数 ， 贪 心 聚 类 算法 
生成 初始 种 群 


<Avg()+Tsve/c >Tavg 


m 


使 用 部 分 映射 交叉 生成 一 个 解 X, 在 调用 
路 径 重 连 算法 寻找 X 与 第 i 个 种 群 的 父 代 蘑 
茵 之 间 更 优 的 解 。 更 新 每 个 种 群 的 父 代 以 

及 全 局 最 优 解 


对 种 群 中 每 个 个 体 进行 更 加 细 化 的 邻 
域 搜索 ， 并 记录 搜索 路 径 。 更 新 每 个 
种 群 的 父 代 以 及 全 局 最 优 解 。 


是 -- 扫 疡 迭代 次 数 w 是 否 小 于 最 大 迭代 次 六 


E 


iiie - prec a 
图 7 HDMRO 算法 流程 图 
Fig.7 The flow chart of HDMRO 


3 ” 算 例 测试 与 结果 分 析 
为 了 验证 混合 蘑菇 算法 在 解决 OLRP-FD 模型 的 有 效 性 ， 


本 文 仿照 文献 [1H] 随 机 生成 两 组 规模 不 同 的 算 例 。 其 中 , 客户 


时 车 辆 就 会 产生 


户 与 仓库 之 间 的 距离 ， 或 
纲 剩 余 货 物 量 不 足以 服 
则 车 辆 需要 返回 仓库 进行 装载 ， 然 而 车 辆 从 当 


前 客户 返回 仓库 的 距离 与 之 前 已 经 行驶 的 距离 之 和 超过 最 大 
距离 约束 ， 那 么 仓库 就 只 需要 重新 派出 一 辆 车 为 当前 客户 服 


iml 


务 , 而 此 时 的 额外 距离 则 为 当 六 


机 模拟 程序 的 


体 步 又 如 下 : 


1 客户 到 仓库 之 间 的 距离 )。 随 


a) 生成 每 个 客户 的 实际 需求 。 X We, (a) 在 区 间 d,d] 


内 产生 随机 数 7， 


个 实数 a 在 区 间 [0 
否则 ， 重 新 随机 生成 了 条 


计算 


> 


模拟 所 有 客户 的 实 
b) 沿 着 HDM 
TEM 


需求 计算 ! 


损失 成 本 。 


际 需求 。 


1]; 如 果 acu. 3H 


相应 的 隶属 度 4 ; 
客户 的 实际 需求 为 r。 
HX， 直到 满足 上 述 关系 。 此 程序 将 


(b) 随 机 生成 一 


RO 生成 的 每 条 配送 路 线 移动 ， 根 据 实际 


服务 “失败 ”而 产生 的 额外 距离 以 及 机 会 


的 需求 4= (dd2,q3) 分 别 从 [1,35]，[36,60]，[61,110]， 三 个 区 

间 内 随机 生成 ， 其 他 相关 数据 见 表 1。 本 文 实验 环境 是 在 

Windows 10 系统 下 的 MATLAB R2017a 实现 混合 蘑菇 算法 的 

编码 以 及 相关 算 例 的 实验 。 相 关 参 数 : 父 代 种 群 大 小 N=40; 
子 代 个 数 M=8; 最 大 迭代 次 数 G-600; WIE c-10. 

表 1 测试 算 例 相关 信 息 

Tab. 1 The relative data of tested instances 


客户 ”候选 仓库 数量 /车辆 仓库 仓库 的 启 车 辆 的 启 车 辆 最 大 距 


数量 ”仓库 车 辆 数 Wu 容 成 本 成 本 离 限 制 

30 5/5 300 1200 100 20 300 
50 5/8 500 2000 100 20 300 
3.1 DPI 的 灵敏 度 分 析 


此 部 分 将 DPI 值 设置 以 步 长 0.1 的 长 度 ， 在 区 间 0-1 内 
变化 ， 并 且 ， 为 了 方便 测试 ， 节 约 计算 时 间 ， 将 4P7 值 固定 
为 1( 保 证 仓库 一 定 能 服务 归属 于 它 的 所 有 客户 )。 测 试 数据 如 
表 2、3 所 示 。 

图 7 和 8 分 别 对 应 表 2 和 3 的 内 容 ， 如 图 


7 和 8 所 示 ， 


:202009.00103v1 
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随 着 DPI 值 的 增加 ， 
路 线 距离 长 度 的 


保持 不 变 ， 机 会 


性 降低 ， 车 辆 的 额外 行驶 距离 曾 
损失 成 本 为 0 是 


刘 


启用 的 车 有 


j 数 增加 ， 
增加 ， 那 么 ， 此 时 车 辆 服务 客户 失败 的 可 能 
ti 会 减少 ， 而 仓库 的 开放 成 本 
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同时 导致 计划 行驶 
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足够 ， 而 额外 行驶 距离 


于 每 个 仓库 的 可 


J 车 辆 数 


未 使 车 辆 数 超过 可 | 
个 算 例 中 的 机 会 损失 成 本 为 零 ， 


1 数量 , 所 以 这 
也 表明 算 例 中 每 个 客户 都 


会 被 服务 。 从 图 


7 和 8 可 观察 到 ， 


随 着 DPI 值 的 变化 , 当 API 


值 一 定时 ，DPI 取 0.6 左右 时 ， 目 标 总 成 本 达到 最 小 值 。 


表 2 30 个 客户 算 例 测 试 结 


Tab. 2 The computational results when the number of customer is 30 


DPI “计划 路 线 距 离 额外 行驶 距离 总 行驶 距离 机 会 损失 成 本 仓库 开放 费用 车 辆 费用 ”总 成 本 计算 时 间 /s 
0.1 410.72 156.20 566.92 0.00 300 62 928.92 134.22 
0.2 419.90 134.22 554.12 0.00 300 86 940.12 140.08 
0.3 428.33 132.70 561.03 0.00 300 100 961.03 130.46 
0.4 439.02 123.74 562.76 0.00 300 100 962.76 127.71 
0.5 440.85 123.11 563.96 0.00 300 100 963.96 132.99 
0.6 449.28 39.91 489.19 0.00 300 122 911.19 128.25 
0.7 459.14 16.69 475.83 0.00 300 144 919.83 129.12 
0.8 462.39 0.99 463.38 0.00 300 160 923.38 130.44 
0.9 468.02 0.07 468.09 0.00 300 160 928.09 132.26 
1 495.48 0.00 495.48 0.00 300 190 985.48 133.42 
表 3 50 客户 算 例 测试 结果 
Tab.3 The computational results when the number of customer is 50 
DPI ”计划 路 线 距 离 额外 行驶 距离 总 行驶 距离 机 会 损失 成 本 仓库 开放 费用 车 辆 费 | 总 成 本 ”计算 时 间 /s 
0.1 534.79 219.5 754.29 0.00 300 62 1116.29 195.22 
0.2 538.04 182.28 720.32 0.00 300 86 1106.32 174.30 
0.3 540.73 111.28 652.01 0.00 300 92 1044.01 228.01 
0.4 542.27 104.56 646.83 0.00 300 100 1046.83 210.21 
0.5 551.04 89.76 640.80 0.00 300 104 1044.80 199.77 
0.6 556.31 52.32 608.63 0.00 300 124 1032.63 223.90 
0.7 559.45 35.29 594.74 0.00 300 150 1044.74 222.50 
0.8 570.43 0.61 571.04 0.00 300 162 1033.04 230.61 
0.9 585.09 0 585.09 0.00 300 180 1065.09 243.60 
1 599.46 0 599.46 0.00 300 184 1083.46 174.60 
mE 3.2 算 例 测试 

NP: 于 研究 中 客户 的 需求 是 模糊 的 ， 所 有 客户 和 仓库 的 位 

站 置 都 是 随机 产生 的 ， 同 时 ， 带 模糊 需求 的 开放 性 选 址 路 径 问 

1o obs: 题 模型 的 研究 在 国内 外 很 少 ， 所 以 无 法 将 表 2、3 所 示 的 

ET HDMRO 结果 与 之 前 研究 的 工作 进行 对 比 。 因 此 ， 为 了 评价 

e] raa HDMRO 的 效率 ， 本 文选 择 将 式 (10) 和 (11) 松 弛 ， 使 其 模糊 需 

400 求 的 限制 条 件 变 为 确定 性 需求 ， 与 松弛 后 的 下 界 进行 比较 ， 

woj 。 同时 , 式 (12) 的 不 等 式 限制 要 变 为 等 式 限制 , 松弛 后 的 公式 如 
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Fig.8 The curves of varied costs changed with DPI 


when the number of customer is 30 
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Fig.9 The curves of varied costs changed with DPI 


when the number of customer is 50 


(30) 和 (31) 所 示 。 
Ea 


ieV jeJ 


DaY;< 


jel 


<Q VkeK (30) 


PZ: Viel 


(31) 


显然 ， 当 模糊 需求 的 左边 界 为 客户 的 确定 性 需求 时 ， 客 


户 的 总 需求 相 较 于 需求 模糊 时 会 下 降 ， 同 时 ， 仓 库 和 车 辆 的 


利用 率 也 会 降低 。 


一 列 为 算 例 规模 ， 
运行 时 间 对 比 ， 其 中 ， 


如 表 4 所 示 ，CPLEX 求解 


Gap = 


因此 ， 可 以 得 到 确定 性 需求 下 的 总 成 本 
相对 应 的 模糊 需求 下 总 成 本 的 一 个 下 界 的 结论 。 表 4 中 ， 
接 下 来 三 列 为 解 的 质量 对 比 ， 最 后 三 列 


是 
第 
为 


的 解 -CP 求解 
HDMRO 的 解 - C LEX 求 解 下 界 .100m 


CPLEX 求 解 下 界 ° 
的 两 个 测试 实例 的 下 界 分 别 是 


503.08 和 677.31， 第 


这 种 使 / 


个 算 例 所 | 
个 算 例 在 3 小 时 限制 内 未 得 到 最 优 答案 。 而 HDMRO 求解 
个 算 例 的 时 间 分 别 为 128.25 秒 和 223.90 $, 
精确 算法 的 求解 器 相 比 ， 该 算法 在 计算 时 间 上 


时 间 为 4011.34 秒 ， 第 二 


y 
因此 , E; CPLEX 


更 高 的 效率 , 可 以 在 较 短 的 计 


算 机 运行 时 间 内 求 得 满意 的 解 。 
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表 4 测试 算 例 结 
Tab.4 The results of tested instances 
算 例 规模 解 的 质量 运行 时 间 /s 
HDMRO(DPI=0.6 和 API-1) CPLEX 求解 下 界 Gap(%) HDMRO CPLEX 12.8 Gap(96) 
30-5 911.19 503.08 81.1296 128.25 4011.34 3027.75 
50-5 1032.63 677.31 52.4696 223.90 10800 4723.58 


除 此 之 外 ， 为 了 进一步 验证 HDMRO 和 OLRP-FD 模型 软件 下 3 小 时 内 求解 出 的 下 界 ， 相 对 于 本 文 使 用 HDMRO K 
的 有 效 性 ， 本 文选 取 了 标准 CLRP 算 例 库 (http://prodhonc. 解 带 模糊 需求 的 OLRP 结果 来 说 ，HDMRO 在 求解 此 类 复杂 
free.fr/Instances/instances us.htm) 中 的 7 组 标准 算 例 进行 相应 问题 上 有 着 较 好 的 效果 ， 除 此 之 外 ，Gap 值 在 7 组 算 例 中 至 
改造 来 进行 实验 。 改 造 步 又 如 下 : 少 大 于 19.48%, 这 说 明了 总 成 本 会 因为 不 确定 的 决策 情况 会 

a) 标准 算 例 的 默认 需求 为 模糊 需求 的 左 界 (d), 右 界 (d3) ”导致 巨大 的 差异 ， 所 以 ， 当 决策 者 所 拥有 的 用 于 决策 的 信息 
为 左 界 的 3 倍 ， 而 da J3[1.5 di, 2.541] 中 随机 生成 一 个 整数 。 不 足 时 ， 利 用 科学 的 手段 进行 合理 的 预测 和 谨慎 的 决策 能 在 

b) 比较 车 辆 容量 2 与 所 有 客户 中 最 大 的 ds. ROCK 一 定 程度 上 避免 企业 损失 (此 组 算 例 中 客户 的 机 会 损失 成 本 
d3， 则 默认 车 辆 容量 不 变 ， 否 则 ， 将 最 大 的 d; 设置 为 车 辆 容 ”为 所 对 应 客户 实际 需求 的 20 倍 )。 


量 。 仓 库容 量 了 设置 为 默认 值 的 三 倍 。 其 余 参 数 使 用 数据 库 为 了 更 进一步 验证 HDMRO 的 求解 效率 ， 本 文选 择 了 
默认 参数 ， 不 做 改变 。 CLRP 标准 数据 库 中 的 16 组 算 例 进行 比较 , 使 用 HDMRO 运 


表 5 中 给 出 了 7 组 改造 后 的 OLRP 实例 在 CPLEX 优化 。 行 每 组 算 例 10 次 后 取 其 中 最 优 值 的 比较 结果 如 表 6 所 示 。 
表 5 HDMRO 在 7 组 标准 OLRP 实例 上 的 计算 结果 
Tab.5 The computational results of HDMRO in 7 standard instances 

HDMRO 的 解 (DPI= 0.6 and API= 1) 


算 例 OLRP DER NNI . Gap(%) 
计划 行驶 距离 额外 行驶 距离 仓库 开放 费 车 辆 费用 总 费用 
20-5-la 43958 20700 6626 25549 12000 64875 47.58 
20-5-lb 34497 17900 5403 15497 5000 43800 26.97 
20-5-2a 41125 19100 6431.5 22769 12000 60300.5 46.63 
20-5-2b 32520 19100 563 13911 6000 39574 21.69 
50-5-la 78562 61400 10825 15385 24000 128772 63.91 
50-5-lb 63099 42000 11680 15385 12000 85454 35.43 
50-5-2a 73242 51600 6324 29319 29000 133065 81.67 


K6 HDMRO 与 其 他 三 种 算法 比较 结果 
Tab.6 The comparative results of HDMRO and other three algorithms 
GRASP 运行 时 间 /s MAPM 运行 时 间 /s LRGTS 运行 时 间 /s ”2-phase HGTS ”运行 时 间 /s HDMRO ”运行 时 间 /s 


Gaskell67-21x5 429.6 0.2 424.9 0.3 424.9 0.2 424.9 6 424.9 13.7 
Gaskell67-22x5 585.1 0.2 611.8 0.3 587.4 0.2 585.1 9 585.1 20.5 
Gaskell67-29x5 515.1 0.4 512.1 0.8 512.1 0.4 512.1 11 512.1 22.3 
Gaskell67-32x5 571.9 0.6 571.9 0.8 584.6 0.6 562.2 40 568.5 19.8 
Gaskell67-32x5 504.3 0.5 534.7 1.0 504.8 0.5 504.3 22 504.3 20.5 
Gaskell67-36x5 460.4 0.8 485.4 1.4 476.5 0.7 460.4 39 479.5 18.3 

Min92-27x5 3062.0 0.4 3062.0 1.0 3065.2 0.3 3062.0 11 3082.0 20.2 
Christofides69-50x5 599.1 23 565.6 3.8 586.4 24 580.4 45 565.6 30.8 
20-5-la 55021 0.2 54793 0.3 55131 0.4 54793 3 54793 15.7 
20-5-1b 39104 0.2 39104 0.3 39104 0.2 39104 4 39104 12.8 
20-5-2a 48908 0.1 48908 04 48908 0.5 48945 3 48908 132 
20-5-2b 37542 0.2 37542 0.3 37542 0.1 37542 4 37542 12.9 
50-5-1 90632 1.8 90160 2.6 90160 0.3 90402 27 90160 30.7 
50-5-1b 64761 1.8 63242 32 63256 1.0 64073 27 66468 352 
50-5-2 88786 24 88298 34 88715 1.8 89342 23 89809 312 
50-5-2b 68042 2.5 67893 2.9 67698 2.0 68479 21 67340 36.7 

取得 最 优 值 次 数 7 0 6 10 11 


表 6 的 16 组 测试 算 例 中 ，GRASPP51( 贪 心 随机 自 适 应 搜 编码， 在 Dell PC Optiplex GX260 装载 2.4 GHz Pentium 4 的 
索 算 法 ) 能 够 求 得 7 组 最 优 解 ，MAPMP9( 基 于 人 群 管理 的 文 ”CPU 和 512MB RAM 以 及 Windows XP 系统 下 进行 测试 ，2- 
化 基因 算法 ) 能 够 求 得 10 组 最 优 解 ,LRGTSP7( 基 于 拉 格 朗 phase HGTS 算法 是 Linux 11.04 下 2GBRAM 以 及 intel core 
松弛 的 禁忌 搜索 算法 ) 能 够 求 得 6 组 最 优 解 ，2-phase ”Duo(2.0GHz) 的 配置 下 用 C++ 编写 的 ), 无 法 进行 直观 的 比较 ， 
HGTSPESI( 两 阶段 启发 式 散 点 禁忌 搜索 ) 能 够 求 得 10 组 最 优 解 ， 但 可 以 肯定 的 是 ，HDMRO 能 够 在 可 接受 的 时 间 内 求解 出 质 
而 HDMRO 能 够 求 得 11 组 最 优 解 ， 相 比 于 另外 四 种 算法 ， 量 令 人 满意 的 解 ， 所 以 ， 根 据 以 上 数据 可 以 得 到 HDMRO 算 
HDMRO 能 够 求解 出 更 多 的 最 优 解 ， 在 求解 质量 上 具有 更 大 ”法 在 此 类 问题 的 求解 有 具有 良好 的 效果 。 
的 优势 。 此 外 ， 在 求解 时 间 这 一 方面 ， 由 于 计算 机 硬件 与 编 ” 3.3 OLRP 与 CLRP 对 比 
译 软件 的 不 同 (文章 中 所 引用 的 前 三 种 算法 都 是 由 C++ 进行 此 部 分 选择 改造 后 的 7 组 算 例 进行 比较 ， 分 析 开 放 式 选 
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址 路 径 问 题 与 普通 选 址 路 径 问 题 以 及 两 者 在 需求 模糊 时 的 区 机 会 损失 成 本 会 让 总 费用 增加 ， 但 OLRP-FD 相对 于 CLRP- 


别 ， 判 断 企 业 是 否 选择 配送 服务 外 包 。 此 部 分 所 有 算 例 使 用 
HDMRO 进行 求解 。 相 关 参 数 : 需求 模糊 时 DPI = 0.6, API 
=1， 车 辆 启用 费用 为 1000， 单 位 距离 费用 为 100， 其 余 参 数 


FD 增加 的 更 少 , 因此 , 企业 就 有 更 多 费用 空间 去 考虑 物流 服 
务 是 否 承包 给 第 三 方 物流 服务 企业 。 其 次 ， 当 客户 规模 为 20 
时 ，Gap2 (EKF Gapl 值 的 原因 ， 主 要 是 因为 此 组 算 例 中 客 


设置 与 3.2 中 的 数据 相同 。 求 解 结果 如 表 7 所 示 。 


CLRP - OLRP 3 " 
— — x100% , 需求 是 JT 
cup xl00% ， 在 需求 是 确定 的 情 


Gapl= ÉF, OLRP 5 


户 较 少 ， 车 辆 服务 所 分 配 的 客户 时 并 未 达到 最 大 行驶 距离 限 
制 ， 能 够 顺利 的 服务 完 所 有 客户 ， 仓 库 无 须 启用 更 多 车 辆 
从 而 没有 产生 机 会 损失 成 本 ， 同 时 ， 在 考虑 额外 行驶 距离 


b ~ 
cr 


CLRP 之 间 的 Gapl 至 少 为 15.44%， 这 说 明 在 不 考虑 最 后 一 ”情况 下 ，CLRP-FD 的 产生 的 额外 距离 高 于 OLRP-FD 产生 的 
段 运输 距离 ， 假 设 企业 的 任何 其 他 与 此 物流 运输 相关 的 费用 额外 距离 , 所 以 Gap2 值 会 大 于 Gapl 值 。 而 在 客户 规模 为 50 
和 第 三 方 物流 公司 的 费用 一 致 的 情形 下 ， 决 策 者 需要 考虑 的 。 时 ,Gap2 值 小 于 Gapl 值 ,是 由 于 最 大 距离 限制 以 及 HDMRO 
是 承包 给 第 三 方 物流 公司 的 服务 费用 与 此 段 距 离 所 节省 费用 ”对 路 线 规划 的 原因 ， 仓 库 启用 的 车 辆 数 有 限 ， 从 而 导致 没有 
的 差 值 的 大 小 ， 由 此 来 决定 是 否 将 物流 配送 服务 外 包 。 而 在 足够 的 车 辆 服务 完 归 属于 该 仓库 的 客户 ， 此 时 会 产生 较 大 的 
需求 是 不 确定 的 情形 下 时 ，OLRP-FD 5 CLRP-FD 之 间 的 。 机 会 损失 成 本 ， 虽 然 在 总 成 本 中 占 主要 部 分 的 还 是 计划 行驶 
CLRP-FD - OLRP-FD . FT 路 线 成 本 , fH OLRP-FD 的 机 会 损失 成 本 高 于 CLRP-FD 的 机 
Gop2( Gap2- OLRPFD “0 )EPA 19.660% 这 说 明 。 会 损失 成 本 导致 Gap 值 有 一 定 的 减 小 ， 从 而 导致 Gap2 值 小 
在 需求 是 不 确定 的 时 候 ，CLRP-FD 的 总 费用 会 比 OLRP-FD — T Gapl fü. 
更 大 ， 那 么 说 明 考 虑 模糊 需求 时 ， 虽 然 车 辆 额外 行驶 距离 与 
表 7 OLRP 与 CLRP 结果 对 比 
Tab.7 The comparative results of OLRP and CLRP 
算 例 OLRP CLRP Gapl(96) OLRP-FD CLRP-FD Gap2(%) 
20-5-1a 43849 54793 24.96 64875 104242 60.68 
20-5-1b 33564 39104 16.51 43800 52413 19.66 
20-5-2a 41125 48908 18.93 60300.5 102347 69.73 
20-5-2b 32520 37542 15.44 39574 52605 32.93 
50-5-1a 63985 90160 40.91 128772 176215 36.84 
50-5-1b 51485 66468 29.10 85454 106676 24.83 
50-5-2a 67519 89809 33.01 133066 176310 32.50 
4 结束 语 D]. 系 d din 27 0%: ipe (Zheng Bin, Ma Zujun, PE 
Fuzzy multi-objective location-routing problem of emergency logistics 
本 文 对 带 模 糊 需 求 的 开放 式 选 址 路 径 问 题 建 立 以 仓库 开 system [J]. System Engineering, 2009, 27 (08): 21-25.) 


放 、 使 用 车 辆 、 车 辆 行驶 距离 、 et 额 
外 行驶 距离 等 成 本 之 和 最 小 为 目标 建立 数学 模型 。 同 时 ， 提 
出 了 求解 离散 问题 的 混合 蘑菇 繁殖 算法 ， 使 用 部 分 映射 交叉 
和 路 径 重 连 算法 组 合 代替 算法 优秀 父 代 蘑菇 细 化 搜索 空间 与 
I 除 种 群 的 操作 ， 合 理 的 避免 信息 在 离散 和 连续 空间 转换 时 
的 丢失 ， 并 且 能 够 更 好 的 探索 解 的 空间 ， 提 高 算法 的 迭代 效 
率 。 除 此 之 外 ， 邻 域 搜 索 时 使 用 概率 法 对 多 种 邻 域 选 择 ， 加 
强 局 部 搜索 能 力 。 最 后 , 通过 对 DPI 的 灵敏 度 分 析 , HDMRO 
与 CPLEX 求解 改造 算 例 结果 之 间 的 比较 ， 以 及 OLRP 与 
CLRP、OLRP-FD 与 CLRP-FD 求解 结果 之 间 的 对 比分 析 , 分 
析 了 企业 选择 第 三 方 物流 企业 服务 的 可 能 性 ， 验 证 了 算法 及 
模型 的 有 效 性 。 
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